Constraints on f{R) theory and Galileons from 
the latest data of galaxy redshift surveys 
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The growth rate of matter density perturbations has been measured from redshift-space distortion 
(RSD) in the galaxy power spectrum. We constrain the model parameter space for representative 
modified gravity models to explain the dark energy problem, by using the recent data of f^{z)as,{z) 
at the redshifts z = 0.06-0.8 measured by WiggleZ, SDSS LRG, BOSS, and 6dFGRS. We first 
test the Hu-Sawicki's f{R) dark energy model, and find that only the parameter region close to 
the standard A Cold Dark Matter (ACDM) model is allowed (A > 12 and 5 for n = 1.5 and 2, 
respectively, at 95% CL). We then investigate the covariant Galileon model and show that the 
parameter space consistent with the background expansion history is excluded by the RSD data at 
more than IOct because of the too large growth rate predicted by the theory. Finally we consider 
the extended Galileon scenario, and we find that, in contrast to the covariant Galileon, there is a 
model parameter space for a tracker solution that is consistent with the RSD data within 2a level. 
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I. INTRODUCTION 

The observational support for the existence of dark 
energy has motivated the idea that the gravita- 

tional law may be modified from General Relativity 
(GR) at large distances to realize the late-time cos- 
mic acceleration (see Refs. 0-01 for reviews). In this 
vein many modified gravity models of dark energy have 
been proposcd-including those based on f{R) gravity 
scalar-tensor theories Q, the Dvali-Gabadadze-Porrati 
(DGP) braneworld scenario Galileons [O, [l^li and 
so on. 

Measuring the growth rate of large scale structures in 
the Universe gives a strong test for the modified gravity 
scenario, because modified gravity models generally pre- 
dict different growth rates from that in the standard A 
Cold Dark Matter (ACDM) model. In fact the equations 
of linear matter density perturbations have been derived 
for a number of modified gravity models-including f{R) 
gravity [lllil, the DGP model [H, Galileons [lij. Re- 
cently De Felice et al. [l3] derived the full perturbation 
equations as well as the effective gravitational coupling to 
non-relativistic matter in the most general scalar-tensor 
theories in 4 dimensions [l8l - l20j (which cover most of the 
single- field dark energy models proposed in literature). 

One of the methods to measure the cosmic growth rate 
is redshift-space distortion (RSD) that appears in clus- 
tering pattern of galaxies in galaxy redshift surveys be- 
cause of radial peculiar velocities. RSD on large and 
linear scales reflects the velocity of inward collapse mo- 
tion of large scale structure, which is directly related to 
the evolutionary speed of matter overdensity perturba- 
tions simply from the mass conservation |21| . Recent 
galaxy redshift surveys have provided measurements of 
the growth rate /m(^) or/m(2)(78(z) as a function of 
redshift up to z ~ 1 |22h3(]( . where fm = d\nSm/dlna, 
dm is the fractional over-density of non-relativistic mat- 



ter, a is the scale factor of the Universe, and as is the 
rms amplitude of over-density at the comoving 8h~^ 
Mpc scale {h is the normalized Hubble parameter Hq = 
100 /i kmscc~^Mpc~^). In this paper we constrain some 
of the modified gravity models proposed to solve the dark 
energ y pr oblem, by using the latest RSD data of the 2dF- 
GRS m, WiggleZ 1^, SDSS LRG [13, BOSS [H, and 
6dFGRS m. 

In this paper we choose the f{R) gravity and Galileon 
models as the representative theories of modified gravity. 
A general difficulty of dark energy models based on modi- 
fied gravity is the emergence of "the fifth force" that can 
violate the constraints from local gravity experiments. 
There are two mechanisms to suppress the propagation 
of the fifth force in local regions: (i) chameleon mecha- 
nism [Slj, and (ii) Vainshtein mechanism [s^. The f{R) 
gravity and Galileons are representative models to avoid 
the fifth force problem by the mechanisms (i) and (ii), 
respectively (see Sec. |lTl where we give a brief review of 
these theories). 

In f{R) gravity several authors put observational 
bounds on the parameter B = (FrR/F){H/H) by 
assuming the ACDM background 33l-[35|. where F = 
df/dR, F^n = dF/dR, and H is the Hubble parame- 
ter (see also Refs. [33 for related works). At the level 
of perturbations the parameter B characterizes the de- 
viation from the ACDM model. The joint data analysis 
of cluster abundance, the cosmic microwave background 
(CMB), and other observations shows that the value of 
B today is constrained to be < 1.1 x 10"'^ at the 
95 % confidence level (CL) [35}. The matter and the ve- 
locity power spectra were also computed with A^-body 
simulations for the Hu-Sawicki model (a popular and vi- 
able f{R) dark energy model, see Eq. ^ below) with 
n = 1/2 ■ However, those past works did not place 
explicit constraints on the viable parameter space of the 
Hu-Sawicki model from the RSD data. 

Recently the covariant Galileon dark energy model 
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(whose cosmological dynamics was studied in Refs. |38l - 
lloj ) was confronted with observations [4l| by using the 
RSD data of WiggleZ and BOSS as well as the data of 
Supernovae la (SN la), CMB, and Baryon Acoustic Os- 
cillations (BAO). It was found that this model is severely 
disfavored over the ACDM. In this paper we show that 
the covariant Galilcon is indeed excluded at more than 
lOcr CL by using the most recent RSD data of SDSS LRG 
and 6dFGRS in addition to the WiggleZ and BOSS data. 
We then further consider the extended Galileon scenario 
in which a tracker solution with an arbitrary constant 
dark energy equation of state smaller than — 1 is realized 
during the matter era. Unlike the covariant Galileon we 
show that there are some viable parameter spaces com- 
patible with the current observational data of RSD as 
well as SN la, CMB, and BAO. 

This paper is organized as follows. In Sec. |IT] we re- 
view the basic properties of dark energy models based on 
f{R) gravity and Gahleons. In Sec.|lll]we study the evo- 
lution of the growth rate of matter perturbations in these 
models. In Sec. IIVI we put observational constraints on 
the parameters of each model, and Sec. |V] is devoted to 
conclusions. 



II. MODIFIED GRAVITY MODELS 

In general the modified gravity models of dark energy 
are required to recover the Newton gravity at short dis- 
tances for the consistency with local gravity experiments 
in the solar system (4^ . As we mentioned in Introduc- 
tion, there are two known mechanisms for the recovery 
of the Newton gravity in local regions. 

One is the chameleon mechanism (3l| , under which the 
mass of a scalar degree of freedom is different depending 
on the matter densities in the surrounding environment. 
If the effective mass is sufficiently large in the regions 
of high density, the coupling between the field and non- 
relativistic matter can be suppressed by having a thin- 
shell inside a spherically symmetric body. 

In f{R) gravity there exists a scalar degree of freedom 
("scalaron" [i^) with a potential of gravitational origin 
for nonlinear functions / in terms of the Ricci scalar R. In 
this case it is possible to design the scalar potential such 
that the chameleon mechanism works at short distances 
by choosing appropriate forms of f{R) [13431 ■ Explicit 
f{R) models of dark energy that can satisfy both local 
gravity and cosmological constraints have been proposed 
in Refs. g^^lll. 

Another mechanism for the recovery of Newton gravity 
is the Vainshtein mechanism (32j . under which nonlinear 
scalar-field self interactions can suppress the propagation 
of the fifth force at short distances even in the absence of 
the field potential. In the DGP model, a brane-bending 
mode (f> gives rise to the self-interaction of the form XD<f) 
(where X = — (V</')f/2) through the mixture with a 
transverse graviton |52l |. This self-interaction leads to 
the decoupling of the field (f> from matter within a radius 



much larger than the solar-system scale 1531. However 
this model is plagued by the ghost problem [5J| , in addi- 
tion to the incompatibility with cosmological constraints 
at the background level |55| . 

The field self-interaction XD<j>, which is crucial for 
the recovery of GR in local regions, gives rise to the 
field equations invariant under the Galilean shift (t>{x) — >■ 
(j){x) + bfj^x^ -|- c in flat spacetime. Nicolis et al. [ll| 
derived the general field Lagrangian by imposing the 
Galilean symmetry in Minkowski spacetime. Deffayet et 
al. obtained the covariant version of the Galileon 

Lagrangian keeping the field equations up to second or- 
der in curved backgrounds (see Eq. pT|) below). In this 
case the Galilean symmetry is recovered in the limit of 
Minkowski spacetime. For the covariant Galileon the 
presence of field self interactions other than X\I\(j) gives 
rise to the theoretically allowed parameter space in which 
ghosts and instabilities are absent [s^ . There exists 
the allowed parameter region in which the background 
cosmological constraints are satisfied . 

If we extend the Galileon self-interaction XD^ to 
the form X^D^ [s^], it gives rise to a scenario equiva- 
lent to the Dvali- Turner model at the background level 
[sst . In this case the dark energy equation of state is 
wde = ^1 ^ l/(2p — 1) during the matter era. How- 
ever there is an anti-correlation between the large-scale 
structure (LSS) and the Integrated-Sachs- Wolfe (ISW) 
effect in CMB, which provides a tight bound on p, as 
p > 4.2 X 10^ (95 %CL) [111, in which case the model is 
practically indistinguishable from the ACDM. This situ- 
ation is alleviated by taking into account other extended 
Galileon terms 0|. In this extended Galileon scenario 
there exist viable parameter spaces in which the LSS and 
ISW are positively correlated [6l|, so that wde docs not 
need to be extremely close to — L 

In what follows we review the dark energy models 
based on f{R) gravity, covariant Galileon, and extended 
Galileon, respectively. 



A. f{R) gravity 

The first class is f(R) gravity [i^, [6^ given by the 
action 

S=^ j d^X^gf{R)+Sra, (1) 

where A/pi is the reduced Planck mass, g is a determi- 
nant of the spacetime metric g^^, f{R) is a function of 
the Ricci scalar i?, and Sm is the action of non-relativistic 
matter. We assume that non-relativistic matter (with a 
negligible pressure) is described by a barotropic perfect 
fiuid minimally coupled to gravity. We focus on the met- 
ric formalism under which the action (jlj is varied with 
respect to g^^ (see Refs. for reviews). 

On the flat Friedmann-Lemaitre-Robertson- Walker 
(FLRW) background with the line-element ds^ — ~dt^ + 
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a^{t)dx'^ the field equations following from the action ([T]) 
are Hi 



by 



3FH^ = {FR- f)/2-3HF 



-2FH = F-HF 



(2) 
(3) 



where F = df/dR, H = d/a is the Hubble param- 
eter, is the energy density of non-relativistic mat- 
ter, R = 6{2H^ + H), and a dot represents a deriva- 
tive with respect to cosmic time t. From these equa- 
tions we can also obtain the continuity equation pm + 
3Hpm ~ 0. Defining the density parameters JIde = 
{FR-f)/{6FH^)-F/{HF) and = p^/i3FH^M^,), 
it follows that fioE + = 1 from Eq. ([2]) • 

When we apply /(i?) gravity to dark energy the model 
is usually constructed to possess dc Sitter solutions char- 
acterized by i? = i?ds = 12i7jg = constant, where H^s 
is the Hubble parameter at the de Sitter solution. Using 
Eq. © with p„i = 0, it follows that 



RdsF{Rds) — 2/(i?ds) : 



(4) 



where i?ds is the similar order to the Hubble parameter 
Hq today. Considering the homogeneous perturbations 
around R = RdS, the de Sitter solution is stable under 
the following condition j63l - l65| 



< RF^b{R)/F{R) < 1 , 



(5) 



in which case the solutions with different initial condi- 
tions converge to the de Sitter attractor. 

For the construction of viable f{R) dark energy models 
we require the condition F{R) > to avoid the appear- 
ance of ghosts [H . The mass squared of a scalar degree 
of freedom is given by ~ (3Fi^(i?))~^ in the regime 
» R, so that the condition > needs to be 

satisfied to be free from tachyonic instabilities [l^. In 
the regions of high density {p HqM^^) the mass M 
needs to be much larger than Hq , so that the chameleon 
mechanism can be at work to suppress the propagation 
of the fifth force. 

A typical example satisfying the above demands is 
given by [13] 



fiR) = R-XRc 



{R/Rc 



\2n 



[R/RcY^-^ + 1 ' 



(6) 



where n, A, Rc are positive constants. There exist other 
models such as f{R) ^ R - \R,[l - (1 + R'^/Rl)~"] 
(n > 0) mi and f{R) = R - \Rctanh{R/Rc) jH^, but 
the basic properties are similar to those in the model © 
(the latter corresponds to the limit n ^ IV The model 
f{R) = R-XRc{R/RcY (0 < p < 1) [^il is also cos- 
mologically viable, but local gravity constraints put the 
severe bound p < 10^^° [1^ (i.e., indistinguishable from 
the ACDM model). For the model (|6]) the modification 
of gravity manifests itself even for n of the order of unity. 

In the following we focus on the model ([5]). The con- 
dition (U) for the existence of de Sitter solutions is given 



7^2"-^(2 + 27^2»-2n)' 



(7) 



where 72.ds = Rds/Rc- The stability condition ([5]) reads 

27^ds-(2n-l)(27l + 4)7^ds + (2n-l)(2n-2) > 0. (8) 

For given n this gives the lower bounds on TZds a-nd A. If 
ri = 1, for example, one has TZds > %/3 and A > 8\/3/9. 
For A and n of the order of unity, Eq. ([7]) shows that i?ds 
is roughly of the order of Rc- In the limit that R ^ Rc 
the model ^ behaves as f{R) ~ R-XRc[l + iR/ Rc)~^% 
which is close to the ACDM model. 

In the regions of high density the chameleon mecha- 
nism works for the model ([ 6l) , in which case the fifth force 
was evaluated in Refs. [45|, From the local gravity 
constraints using the violation of weak equivalence prin- 
ciple, we obtain the following bound [45l | 



n > 0.9. 



(9) 



Note that n and A arc the two fundamental parame- 
ters for the model ([5]). For given those parameters the 
ratio TZds = Rds/Rc is known from Eq. ([7]), which leads 
to the estimation Rc « 11^/71^8 ■ Theoretically natural 
model parameters are n and A of the order of unity, but 
it remains to see whether such parameter space can be 
consistent with the RSD data. 



B. Covariant Galileon 

The second class of dark energy models is the covari- 
ant Galileon [ll|, [T^j . This theory is the sub-class of the 
Horndeski's general action p^l - l20l | characterized by 
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(10) 



where 



C2^K{<p,x), /:3 = -G3(0,A)n<^, 

£4 = G4(0, X)R + G4.X [(□</>)' - (V^V,</)) (V'V'' </))], 

£5 = G5(0,A)G,,, (V^V^^) 

-(Gs.jf /6) [{n<j,f - 3(00) (V^V,0) (V^V^) 
+2(V'^V„0) (V"V^0) (V^V^^)] . (11) 

Here K and Gi {i — 3, 4, 5) are functions in terms of a 
scalar field cj) and its kinetic energy X = —g^'^ d^(j)du(t>/2, 
and G^i, is the Einstein tensor. The covariant Galileon 
[T^ corresponds to the choice 

K = ciM^(t)-C2X, 



G3 = CgA/M^ 

^pl 

G5 = ic^X^/M'^ 



Gi = MlJ2 - caX^/M^ , 



(12) 
(13) 
(14) 

(15) 
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where c^'s are dimensionless constants, and M is a con- 
stant having a dimension of niass^. Since we are inter- 
ested in the case where the cosmic acceleration is driven 
by field kinetic terms without a potential, we set ci = in 
the following discussion. Note that the Einstein-Hilbert 
term appears inside the Lagrangian £4. 

On the flat FLRW background the equations of motion 
for the covariant Galilcon are [s^ 

SM^^H^ = puE + Pm , (16) 
3Mp2 JJ^ + 2M^iH = -Pde , (17) 

where 

PUE = -C20V2 + 3c3-ff0VM3-45c4-ff20V(2Af^) 
+21c5H^(j>yM^, (18) 

Pde = -C2</>V2 - C3(A'<^/m3 

+3c4<p^[8H4> + (3ff2 + 2H)^]/{2M^) 
-3c5H^'^[5H4> + 2{H^ + H)<})]/M^ . (19) 

The dark energy equation of state wde is defined as 
wde = -Pde /pde- 

The de Sitter solution {H = H^s = constant) is real- 
ized for 4> = (pds = constant. Normalizing the mass M 
to be = J^IpiHls and defining a;ds = ipds / {HdsMpi) , 
Eqs. ([TB)) and P7)) lead to the following relations at the 
de Sitter solution [s^ 

C2xls = 3(2 + 3a -4/3), (20) 

C3xls = 2 + 9(a-/3), (21) 

where 

a = C4x'^s ' ^ = csa^ds • (22) 

The relations ((20)l and ((2T|) are not subject to change 
under the rescaling Xds — > JXds, C2 —J- 22/7^, and 
C3 — > C3/7'', where 7 is a real number. Hence the rescaled 
choices of C2 and C3 lead to the same cosmological dynam- 
ics. This means that we can set Xds = 1 without loss of 
generality, in which case C2 and C3 are directly known 
from Eqs. ([^0]) and (HI]). It is convenient to use the two 
model parameters a and f3 because the coefficients of 
physical quantities can be expressed in terms of them. 

Introducing the two variables ri = (f>dsHds/i4'H) and 
^2 = {4'/4'ds)'^ /ti, which satisfy ri = r2 = 1 at the de 
Sitter solution, the dark energy density parameter fioE = 
PDE/(3Mp[iJ^) can be expressed as 

f^DE = -{2 + 3a-4l3)rlr2/2+{2 + 9a~9l3)rlr2 

-15arir2/2 + 7pr2. (23) 



^ In the original paper of Galileons [Till the linear potential ciM^ip 
is expressed as a separate Lagrangian Ci . Here we include both 
the potential and the kinetic term inside £2- 



The autonomous equations for the variables ri and r2 
were derived in Ref. [401. There are three distinct fixed 
points for this dynamical system: (A) (ri,r2) = (0,0), 
(B) (ri,r2) = (1,0), and (C) {n,r2) = (1,1). The 
point (C) is the de Sitter solution, which is always stable 
against homogeneous perturbations. The point (B) cor- 
responds to a tracker solution characterized by cx 1/ H. 
For the initial conditions close to the fixed point (A) the 
solutions converge to the tracker once ri approaches 1. 
The fixed point (B) is followed by the de Sitter point (C) 
after r2 grows to the order of 1. The epoch at which 
the solutions approach the tracker depends on the initial 
values of ri. For smaller ri the tracking occurs later. 

Along the tracker one has f^DE = ?'2 and wde — 
—2/(1 -I- JIde) after the radiation-dominated epoch. In 
the regime close to the fixed point (A) the dark energy 
equation of state is given by wue = ^1/8 during the 
matter era. After the solutions reach the tracker, wbe 
changes from —2 (matter era) to —1 (de Sitter era). 

The joint data analysis of SN la, CMB, and BAO shows 
that the tracker is disfavored from the observational data 
of the distance measurement (56| . For the compatibility 
with the data the solutions need to approach the tracker 
at late times with the minimum value of wde larger than 
— 1.3. Taking into account the conditions for the avoid- 
ance of ghosts and Laplacian instabilities j39l . lioj . the 
parameters a and /3 are constrained to be 

a 1.404 ± 0.057, /? = 0.419 ± 0.023 , (24) 

from the joint data analysis of SN la (Union2), CMB, 
and BAO 

We stress that, for given initial conditions of ri and 
r2 (in other words, given initial conditions of (j) and H) 
under which the solutions enter the tracker at late times, 
there are two parameters a and f3 left for the likelihood 
analysis in Sec.|lVl The bounds ([M)) will be used for such 
an analysis. 

C. Extended Galileon 

The extended Galileon [g^] is described by the Horn- 
deski's Lagrangian pT|) with the choice 

K = -~C2M^^^-P'^ XP' , (25) 

G3 = csM^-^P'XP^ , (26) 

G4 = MpV2 - C4MI^^~^'"'> XP^ , (27) 

G5 = Sc^M'^^+^P^^^XP^^ , (28) 

where c; and pi (i = 2, • • • , 5) are dimensionless con- 
stants, and A'/i's are constants having a dimension of 
mass. The covariant Galileon corresponds to p2 = Pa = 1 
and p4 = P5 = 2. In Ref. [60t it was shown that there 
exists a tracker solution with = constant {q > 0) 

for the powers 

P2=P, P3=p+{2q- l)/2 , 

P4=p + 2q, p5=P+(6g-l)/2. (29) 
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Note that the tracker corresponds to an attractor where 
the solutions with different initial conditions converge to 
a common trajectory {p = 1 and q = 1/2 for the covariant 
Galileon) . 

The extended Galileon possesses de Sitter solutions 
with H = Hds ~ constant and = (f>ds = constant. 
We relate the masses Mi's in Eqs. with Has, 

as M, = (i^dsMpl)l/^ M3 = iH~i'''^Ml~'P^y/<^'-^P'^\ 
Mi = (iJ-g2P*M^-2?'^)i/(2-4p4)^ and M5 - 
^^2+2p5^2p5-i)i/(i+4p5). For the existence of de 
Sitter solutions the following relations hold between the 
coefficients 



2 



C2 = o ( — ) (2 + 3a -4/3), 



(30) 



C3 



72 /_2 

2p + 2q-l U'dS 



p+q 



^ [p + 3(p + g)(a-/3)],(31) 



where Xds = 0ds/(^^dsA-^pi), and 



p+2q 



C4 



/3 = 2V2(^, + 3,-i) (^^k^'^'' 



C5 , 



(32) 
(33) 



The background equations on the flat FLRW back- 
ground are given by Eqs. and P7| . with dif- 
ferent forms of pDE and Pde [sO, EH- Introducing 
the dimensionless variables ri — {4>dsl4'Y''HdslH and 

r2 = [(0/0ds)''/?-i]*^+^''^^*^+^''\ the tracker solution cor- 
responds to ri = 1 with the dark energy density param- 
eter f^DE = /0DE/(3i?^A'/pi) = r2. After the radiation era 
the dark energy equation of state wde = -Pde / Pde along 
the tracker is given by [iO, [fiH 



1 



Wde 



1 + sVl 



where s 



DE 



f- (34) 
2q 



During the matter-dominated epoch (fJDE ^ 1), it fol- 
lows that i^DE — — 1 — s. The covariant Galileon corre- 
sponds to s = 1 and wde — —2, in which case the tracker 
is incompatible with the observational data at the back- 
ground level [5^. In the extended Galileon model it is 
possible to realize wde close to —1 for s ^ 1. The joint 
data analysis of SN la, CMB, and BAG shows that the 
parameter s is constrained to be [6]| 



0.034: 



-0.327 
-0.034 



(95 % GL) . 



(35) 



Hence the tracker solution in the range —1.36 < iude < 
— 1 can be consistent with the background cosmology. 

We choose the initial conditions ri = 1 and r2 = 
f^DE 'C 1 in the early matter era. We identify the present 
epoch to be = 0.72, which provides the initial value 
of r2 at a given redshift. For given values of p and q 
there are two parameters a and P left for the likelihood 
analysis in Sec. IIVI 



III. COSMOLOGICAL PERTURBATIONS 
A. Linear Perturbation Equation 

In Ref. [13 the full linear perturbation equations were 
derived for the most general scalar-tensor theories with 
second-order field equations [H, [l^ . In this section we 
study the evolution of the cosmic growth rate in f{R) 
gravity, covariant Galileon, and extended Galileon. 

We consider the scalar metric perturbations ^E* and $ in 
the longitudinal gauge about the flat FLRW background. 
The perturbed line element is then given by [G^I 



ds^ = -(1 + 2^-) dt^ + a^{t){l + 2$)da;2 



(36) 



The energy density of non-rclativistic matter is decom- 
posed into the background and inhomogcneous parts as 
Pm{t) + 5pm(t,x). We write the four velocity of non- 
relativistic matter in the form = (1 — V'w), where 
V is the rotational-free velocity potential. We also intro- 
duce the following quantities 



<5 = 5pm/Pn 



(37) 



In Fourier space with the comoving wavenumber k the 
matter perturbation obeys the following equations of mo- 
tion 



j + 6'/a-h3$ = 0, 

e + HO- (/cVa)* = 0, 



(38) 
(39) 



where k = |fe|. We introduce the gauge-invariant density 
contrast 



5m = S+{3aH/k^)e. 
From Eqs. ^ and ^ it follows that 



2HSm + (fcVa^)* = 3(/ + 2Hi) 



(40) 



(41) 



where / = {aH/k'^)e - $. 

In order to confront the models with galaxy clustering 
surveys, we are interested in the modes deep inside the 
Hubble radius. In this case we can employ the quasi- 
static approximation on sub-horizon scales, under which 
the dominant contributions to the perturbation equations 
are those including the terms fc^/a^ i5, and the mass 
M of a scalar degree of freedom [50, [Sa]- Under this 
approximation we obtain the modified Poisson equation 



(42) 



where GeS is the effective gravitational coupling whose 
explicit form is different depending on the theories [l7| . 
In the framework of General Relativity, G^g is equivalent 
to the gravitational constant G = 1/{8ttM^i). 

Under the quasi-static approximation on sub-horizon 
scales the r.h.s. of Eq. (|4ip can be neglected relative to 
the l.h.s. of it. Since 6^ ~ i5, the matter perturbation 
obeys the following equation 



2H5jn - 4:T:GcSPm.S„ 



0. 



(43) 
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From this equation it is possible to predict 5m, fm, and 
their evolution, and then it can be compared with the 
observed data of fm{z)(Jg,{z). 



B. Comparison with Observations 



The perturbation 5g of galaxies is related to 5m via the 



survey 



bias factor fe, i.e., 5q 



The galaxy power spectrum 



V^ik) in the redshift space can be modelled as j2ll. \24 



Kik) - Vgg{k) + 2y.^Vge{k) + ii^Veeik) 



(44) 



where ^ = k ■ r/{kr) is the cosine of the angle of the k 
vector to the line of sight (vector r). Vggik) and Vgg{k) 
are the real space power spectra of galaxies and 6, respec- 
tively, and Vge{k) is the cross power spectrum of galaxy- 
6 fluctuations in real space, where 6 is related to 6 but 
normalized by the Hubble velocity, i.e., 9 = 9/{aH). In 
Eq. (|44p we did not take into account the non-linear ef- 
fect coming from the velocity distribution of galaxies in 
collapsed structures. 

For the linearly evolving perturbations Eq. (p8| shows 
that 6 is related to the growth rate of matter perturba- 
tions, as 

~e = e/{aH) ~ -fm5m = -frnSg/b , (45) 
fm = 5m/{H5m), (46) 

where we neglected the $ term. In this limit we obtain 
the famous Kaiser's formula 1211: 



Vt{k)^{l+tl'PfVgg{k) 



(47) 



where = fm/b is an anisotropy parameter. There- 
fore, the anisotropy of V'g induced by RSD primarily 
constrains /3, and an independent measurement of bias 
gives a measurement of fm- Instead, (feas)^ can be mea- 
sured from the amplitude Vg , and hence the combination 

P ^gY^'^ gives an estimate of fmC^a that is free from the 
uncertainty of the bias and can be measured only from 
galaxy power spectrum. 

For a given model we can numerically solve Eq. (j43p to 
find 5m.{z) and (5,„(z) at a given scale of x Mpc. In 
GR the cosmic growth rate is independent of scales for 
linear perturbations, but in some modified gravity mod- 
els like f{R) gravity the scale dependence is present even 
at linear level. For the models close to the ACDM in the 
early cosmological epoch, the evolution of matter pertur- 
bations during the deep matter era is given by (5m oc a, 
which corresponds to 5m — H5m and hence /„, = 1. The 
modification of gravity manifests itself in the late cosmo- 
logical epoch. 

Now it is easy to predict fmiz)<7x(z) once the nor- 
malization of ax is fixed, where ax oc 5m is the rms 
fluctuation amplitude at this scale. This is the quan- 
tity that we want to compare with the measurements of 
fm{z)as{z). The quantity fm^s measured from a galaxy 



0.067 
0.17 
0.22 
0.25 
0.37 
0.41 
0.57 
0.6 
0.78 



0.423±0.055 
0.51±0.06 
0.42±0.07 
0.3512±0.0583 
0.4602±0.0378 
0.45±0.04 
0.415±0.034 
0.43±0.04 
0.38±0.04 



6dFGRS (2012) [29] 
2dFGRS (2004) [23] 
WiggleZ (2011) [26] 
SDSS LRG (2011) [27] 
SDSS LRG (2011) [27] 
WiggleZ (2011) [26] 
BOSS CMASS (2012) [28 
WiggleZ (2011) [26] 
WiggleZ (2011) [26] 



Table I: Data of /mcrg measured from RSD with the survey 
references. 



redshift survey is based on (3 and Vg , and hence it should 
essentially be regarded as fm<^x with x being the char- 
acteristic scale of the survey. The conversion into the 
conventional parameter erg is done by assuming the stan- 
dard shape of the linear matter power spectrum. There- 
fore, provided that a given modifled gravity model pre- 
dicts a matter power spectrum similar to that of ACDM 
at z = 0, we can compare the theoretically predicted 
fm{z)ax{z) with the observed fm{z)a^{z), if we choose 
X as the characteristic galaxy survey scale and normalize 
ax{z) as ax{z = Q) = a%{z = Q). For this normalization, 
we use the WMAP7 year bound [toI 



a8(z = 0) = 0.81liro'3?. 



(48) 



The data of /mCg used to constrain the models are 
summarized in Table IH There is another constraint 
/mC[s = 0.49 ± 0.18 at z = 0.77 derived by Guzzo et 
al. 123 , but we exclude this data because of more recent 
Wig gleZ measurement at the redshift z = 0.78 [l^ (see 
Ref . [30| ) . The characteristic scales relevant to these sur- 
veys are 10 /i"^ Mpc < k^'^ < 60 h^^ Mpc. It should be 
noted that, even though we compare the theory and data 
in terms of /mCrg, we are essentially comparing fm<^x- 

Below we discuss the evolution of perturbations in 
three modified gravity models separately, and compare 
them with the observed data for some representative pa- 
rameter sets. 



C. /(;?,) gravity 

The effective gravitational coupling in metric f{R) 
gravity is given by [50| 



Gcff 



G\ + 4,m{k/af/R 
F l + 3m{k/ay/R 



(49) 



where the parameter m = RF^r/ F characterizes the de- 
viation from the ACDM model {f ^ R~ 2 A) . In order to 
avoid ghosts and tachyonic instabilities we require F > 
and Ffl > (with i? > 0), respectively [4^, so that 
TO > 0. 

Equation (P^]) shows that the transition from the "Gen- 
eral Relativistic regime" {GcS = G/F) to the "modified 
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Figure 1: Evolution of fmcrs versus the redshift z for the f{R) model ((6| with n = 1. The left, middle, and right panels 
correspond to the scales = Wh~^, 30 h~^, 60 Mpc, respectively. The solid, long dashed, dashed curves correspond to 
A — 1.55, 10, 20, respectively, whereas the dotted curve corresponds to the ACDM. We also plot the observational bounds on 
Uas constrained from 2dFGRS [H, SDSS-LRG [H, WiggleZ 6dFGRS [H, and BOSS [H. 
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Figure 2; The same as Fig. [T] but for n = 2 and A — 1.55, 5, and 10. 



gravitational regime" (Gofi = iG/{3F)) occurs around 
m ~ {aH/k)^. If such a transition occurs at the redshift 
Zc larger than 1, the transition redshift for the model © 
can be estimated as 1711 



k 



aoHo 



2n(2n- 



1) (2f](?^)2"+i 



\2n 



(0)2("+l) 



l/(6ri+4) 



(50) 

where oq is the today's scale factor normalized as = 1, 
r^Qg and rim' are the today's density parameters of dark 
energy and non-relativistic matter respectively. We iden- 
tify the present epoch to be fi^g = 0.72. The modes with 
larger k enter the modified gravitational regime earlier, 
so that the growth of perturbations tends to be more sig- 
nificant. For smaller n and A the transition redshift also 
gets larger. Recall that n and A are bounded from below 
from Eqs. 0-®. 

Numerically we solve the perturbation equations with- 
out using the quasi-static approximation on sub-horizon 



scales. The initial conditions of the perturbation SR is 
chosen such that the oscillating mode is sub-dominant 
relative to the matter-induced mode (along the lines of 
Rcfs. OHO])- The quasi-static approximation based on 
Eq. ((43|) with the gravitational coupling (|49)) is accurate 
enough to estimate the growth rate of perturbations for 
sub- horizon modes [t^ . [73j . 

In Fig. [T] the evolution of fm<^s is plotted for n = 1 
and three different values of A, with the three different 
wavenumbers of k~^ = 10, 30, and 60 Mpc within 
the characteristic range of galaxy survey scales. When 
A = 1.55 the numerical value of the transition redshift 
is larger than a few and hence the effective gravitational 
coupling in the regime z < 1 is given by Gcff ~ 4G/(3F). 
As we sec in Fig.[Tl /mCs is larger than that in the ACDM 
model. When n = 1 and A = 1.55, the growth rate is 
almost independent of scales, which reflects the fact that 
Zc is larger than 1 on those scales. 

In Fig. [T] we find that the dispersion of /mfs with re- 
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spcct to A appears on larger scales. This comes from the 
fact that, for smaller k and larger A, the transition red- 
shift (|50|) gets smaller. Even for the scale 60 Mpc and 
A = 20, however, /,„(T8 is larger than that in the ACDM, 
so that the f{R) model with n = 1 can be distinguished 
from the ACDM. 

For larger n the transition redshift Zc tends to be 
smaller, and the growth rate of 6m today gets smaller 
for given k and A. In Fig. [2] we show the evolution of 
fm<^8 for 71 = 2 and A = 1.55,5, 10. In this case the dis- 
persion of fmCs with respect to A appears even for the 
scale 10 Mpc. On the scales larger than 30 h'^ Mpc 
the values of /mCs are almost degenerate for A > 10, 
whereas for A = 0(1) the growth rate is still large even 
on the scale 60 h^^ Mpc. 

When compared with the observed data in Figs. [T] and 
[21 the f{R) model ^ with n and A of the order of 1 
is in tension with most of the data. In Sec. IIVI we put 
observational constraints on the model parameters n and 
A. 



0.8 



0.7 



Covariant Galileon/i) «-i-347, M.442 

.1 .1 (ii) a=1. 360, (3=0.433 

k =30h Mpc '\(iii) a=1. 347, (3=0.424 

,-' /^~~^\ ACDM 

/ 2dFGRS 1— •— 



/ / 
/ 

/ / 
• / 
/ 



\ SDSS-LRG 1— e- 
WiggleZ i-^ 




D. Covariant Galileon 



The effective gravitational coupling for the covariant 
Galileon is given by ^ (see also Refs. [TUllIli) 



Goff = 



2Mp2(C|-C3C5) 



-G, 



2C1C2C4 — C1C5 — C2C3 
where the time-dependent coefficients Ci's are 



(51) 



Ci = 2Af 2j + 3c4</,4/Af6 „ 6c5iJ0VA/^ (52) 

C3 = 2M^i~ C4,<p^/M^ ^6c5(j)^^/M^, (54) 
C4 = 12c44>'^4>/M^ + AciH4:^/M^ -6c5H<p^/M^ 



(55) 



C5 = C2 - 4c3H^/M^ ~ 2c30/Af3 + 26c4H^<j)^/M^ 

-36c5H^<j)^4>/M^ - 24:C5H(j)^{H^ + H)/M^. (56) 

At the dc Sitter solution (ri = r2 = 1) the formula 
([5T|) reduces to 



Goff 



G 3(a - 2/3) 



(57) 



This shows that a — 2(3 > to avoid ghosts. The late- 
time tracker corresponds to the solution along which the 
conditions ri <^ 1 and ^2 ^ 1 are satisfied during most 
of the cosmological epoch prior to the dominance of dark 
energy. In this regime, expansion of Gcff around ri = 
and r2 = gives 



~G~ 



1 



255 



13- 



211 



■ari r2 



(58) 



Figure 3: Evolution of /mcrs versus the redshift z for the 
scale 30 Mpc in the covariant Galileon model. The solid, 
long dashed, and dashed curves correspond to the cases (i) 
a = 1.347, /3 = 0.442, (ii) a = 1.360, /3 = 0.433, and (iii) a = 
1.347, P = 0.424, respectively, whereas the dotted curve shows 
the evolution of /mcrg in the ACDM model. The observational 
data are the same as those plotted in Fig. [1] 



The non-ghost condition corresponds to /3r2 > [39[, so 
that GefT > G. Hence the growth rate of 5m is larger 
than that in the ACDM. 

Numerically we solve the full perturbation equations 
presented in Refs. [l^, [l3| without using the quasi- 
static approximation for the model parameters given in 
Eq. (j24|) . The quantity ri needs to be much smaller than 
1 during the early matter era, in which case the solutions 
approach the tracker at late times. The initial conditions 
for the perturbations are chosen to satisfy $ = and 
= in the deep matter era (see Ref. [l^ for detail). 

In order to avoid ghosts and Laplacian instabilities the 
parameters a and 5 are restricted to be in the range < 
a — 2/3 < 2/3 [H, li^l. Numerically we find that, for a — 
2(3 > 0.505, the perturbations show violent instabilities 
during the transition from the matter era to the de Sitter 
epoch. This is associated with the fact that the effective 
gravitational coupling (|5ip tends to cross for larger a. 
Hence we focus on the parameter space < a — 2(3 < 
0.505, in addition to the constraints from the background 
evolution [Eq. (^1)) ]. 

In Fig. [3] we plot the evolution of fmO's for the scale 
k^^ = 30 /i^^ Mpc with several different values of a and 
(3. For larger values of a — 2/3 the growth of Sm tends 
to be more significant, so that fmO's gets larger. Even in 
the case (i) (see caption of Fig. [3]) that gives the smallest 
fmO's within the allowed ranges of a and (3, the model is 
in tension with the recent RSD data because of the large 
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Figure 4: The grey region shows the theoretically allowed pa- 
rameter space for the extended Galileon model with p = 1 
and q — 5/2 (the region in which ghosts and Laplacian 
instabilities are absent). Each border corresponds to (a) 
P = a/2, (b ) 13 = a/2 - 1/15, (c) /3 = (408a + 68 - 
2717^3(272 

4V3630 - 495a) /99, respectively. 



75a)a + 68)/561, and (d) /3 = (242 - 15a + 
Taken from Ref. [sl. 



growth rate of Sm ■ Numerically we also confirm that the 
evolution of fmcrs is practically independent of the scales 
in the range 10 Mpc < k'^ < 60 Mpc. 



E. Extended Galileon 

The effective gravitational coupling of for the extended 
Galileon is given by the form ([ST]) , with more complicated 
coefficients Ci (see Ref. [6l| for detail). At the dc Sitter 
solution (ri = r2 = 1) it follows that 



Goff 



2(1 -p) + 3(l + 2g)(a-2/3) 



(59) 



The general expression of Goff is quite involved even 
along the tracker (ri = 1). In the following we focus 
on the tracker solution (ri = 1) for the model p = 1 and 
q = 5/2, in which case s = 0.2 and hence t«DE = —1.2 
during the matter era. This model shows some deviation 
from the ACDM model, but it is still consistent with 
the background-level observations (sec Eq. (|35|) ). In the 
regime r2 ^ 1, Goff is approximately given by 



Goff 
G 



1 + [48411a2 _ 3q,(3560 + 60291/3) 

+22(50 -I- 924/3 -I- 76Alf3^)]r2 
/[10(308- 1536a -f-3201;3)] , 



whereas at the de Sitter solution Goff / G 



(60) 
l/[9(«-2/3)]. 
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Figure 5: Evolution of /mcrg versus the redshift z for the 
scale 30 h^^ Mpc in the extended Galileon model. The solid, 
long dashed, and dashed curves correspond to the cases (i) 
a = 0.1, /3 = 0.049, (ii) a = 3, /3 = 1.434, and (iii) a = 5, 
P = 2.498, respectively, whereas the dotted curve corresponds 
to the ACDM. The observational data are the same as those 
plotted in Fig. [T] 



In Refs. [60|, |6l| the authors clarified the parameter 
space in which ghosts and Laplacian instabilities of scalar 
and tensor perturbations are absent. Figure U shows such 
a viable region, which is surrounded by the 4 borders 
given in the caption. On the line (a) /3 — a/2, one has 
Goff/G — s- OD at the de Sitter solution and Goff/G ~ 
1 -t- 4(275 - 129a)r2/[5(616 + 129a)] for rs < 1. In the 
latter regime Goff gets larger for smaller a. Hence today's 
values of on the line (a) are large around the region 
a = f3 = 0. The case (i) in Fig. [5] corresponds to such an 
example. 

On the line (b) shown in Fig. 01 we have Goh /G = 5/6 
at the dc Sitter solution and Goff/G ~ 1 + 137*2/25 for 
r2 <C 1. In the regime r2 <C 1, Goff on the line (b) is 
larger than that on the line (a) for the same value of 
a, but at the de Sitter solution Goff is even smaller than 
that for the ACDM. The case (ii) in Fig.[5]corresponds to 
such an example, which exhibits rapid decrease of fmO's 
for z < 0.5. This case shows a better compatibility with 
the RSD data in the low redshift regime relative to the 
case (i). 

Around the line (a) , the effective gravitational coupling 
Goff in the regime r2 <SC 1 is smaller than G for a > 
275/129. Although GgS goes to infinity at the de Sitter 
solution, the growth rates of matter perturbations are 
suppressed in high redshifts for a > 275/129. The case 
(iii) in Fig. [S] corresponds to such an example, which are 
consistent with some of the RSD data because of the 
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suppression of fmO's in the early cosmological epoch. 

We also confirm that the evolution of fmO's is practi- 
cally independent of the scales in the range 10 Mpc 



that there is a severe fine tuning for the initial conditions 
of perturbations to avoid the dominance of the oscillating 
mode [3 [13l • Hence it is not so meaningful to consider 
such cases even if may be reduced relative to the case 
n < 2.5. 



IV. OBSERVATIONAL CONSTRAINTS FROM 
RSD 



B. Covariant Galileon 



Here we determine quantitative model parameter space 
that is excluded from the RSD data for the three modi- 
fied gravity models discussed in Sec. IIIII The confidence 
regions of excluded parameter spaces are calculated by 
the standard analysis. For a given set of model pa- 
rameters, Xobs calculated from the 9 data points of 
fmiz)as{z), and the exclusion confidence level at this 
point is the probability of < Xobs distri- 
bution with d degrees of freedom, where d = Nd — Nm, 
Nd = 9 is the data number, and Nm is the number of 
model parameters. For all the three modified gravity 
models considered here, ~ 2. 



A. f{R) model 

Wc first place observational constraints on the f{R) 
model ([ni for three different scales 10 h~^, 30 h~^, 
60 Mpc. We identify the present epoch to be = 
0.72. In Fig. iniwe plot the parameter space in the (n. A) 
plane constrained by the RSD data shown in Table H] 
The grey regions are forbidden by the violations of the 
de Sitter stability condition (|5]) as well as of the local 
gravity constraint ^ . Figure IH] shows that observational 
constraints tend to be tighter for smaller scales. 

The model with n = 1 and A < 20 is outside the 
95.4 % CL contour for the scales lO/i-^Mpc < k'^ < 
60 Mpc. This reflects the fact that the cosmic growth 
rates in this parameter space are too large to be com- 
patible with the recent RSD data in the regime z < 1, 
irrespective of the values of A (see Fig. [1]). 

For n ~ 2, the parameter A is constrained to be 
A > 8 (lO/i-^Mpc), A > 5 (30/i-iMpc), and A > 3 
(60 h~^Mpc) at 95.4 % CL. In terms of the dimensionless 
parameter m = RF,r(R) / F{R) characterizing the devia- 
tion from the ACDM [6|, we have m{z = 0) = 9.6 x 10"^ 
for n = 2, A = 8 and m{z = 0) = 6.3 x 10"'' for 
n = 2, A = 5. Hence, for n = 2, the today's values of m 
are bounded to be m(z 0) < 9.6 x 10"^ (10/i~iMpc) 
and m{z = 0) < 6.3 x 10"^ (30/i-iMpc) at 95.4 % 
CL. The constraint on m{z — 0) coming from the in- 
tegrated Sachs- Wolfe effect of the CMB anisotropics is 
much weaker, ?7i(z = 0) < 1 [1, H^]- This shows that 
the measurement of the cosmic growth rate from RSD is 
sufficiently powerful to constrain the f{R) gravity model 
tightly. 

If wc choose the values of n greater than 3, the model 
is close to the ACDM. In this case, however, the scalaron 
mass is very heavy in the early cosmological epoch, so 



For the covariant Galileon the allowed parameter space 
is constrained as Eq. from the combined data anal- 
ysis of SN la (Union2), CMB, and BAG. We iden- 

0.73. If the 



(0) 



tify the present epoch to be 51^^ 
Constitution SN la data are used instead of Union2, 
the observational bounds on a and /3 become a = 
1.411 ± 0.056 and ^ = 0.422 ± 0.022 [s^. Set- 
ting the prior a — 2(3 < 0.505 to avoid the neg- 
ative values of Goft, the allowed parameter space is 
the triangle surrounded by the three points (a, /?) = 
(1.347, 0.421), (1.347, 0.442), (1.389,0.442) for Union2 
and {a, (3) = (1.355, 0.425), (1.355, 0.444), (1.393, 0.444) 
for Constitution. In these parameter regions we further 
place constraints on the model from the RSD data. 

In Fig. [7] we show observational bounds on the co- 
variant Galileon model in the (a, (3) plane for the scale 
30 Mpc. The parameter space constrained from the 
background cosmology is excluded at more than llcr CL. 
This comes from the large cosmic growth rate relative to 
that of the ACDM, as shown in Fig. [31 As the model pa- 
rameters approach the border a— 2/3 = 0.505, it is further 
difficult to match with the data because of the significant 
enhancement of matter perturbations. We also find that 
the observational constraints are practically independent 
of the scales. For the scale 60 Mpc the covariant 
Galileon is excluded at more than lOcr CL. 

A similar conclusion was recently reached in the pa- 
per [4l| , in which the authors included the WiggleZ and 
BOSS data as well as SN la, CMB, BAG data. The 
difference from their work is that wc have taken into ac- 
count the 2dFGRS, 6dFGRS, and SDSS-LRG data as 
well. Even though the 2dFGRS and 6dFGRS data pro- 
vide larger values of /mCg relative to those from WiggleZ 
and BGSS, the covariant Galileon cannot be consistent 
with most of the currently available RSD data. 



C. Extended Galileon 

In the extended Galileon model with p = I and q — 
5/2; the theoretically allowed parameter space is given by 
the region shown in Fig. As mentioned in the previous 
section, this plot assumes that the dark energy equation 
of state during the matter era is wde = ^1 ^ s = —1.2, 
where s = p/{2q). This model is interesting because it 
is different from the ACDM model but is still consistent 
with the joint constraint of SN la, CMB, and BAG [6l| 
(see also Eq. (pSjl ). The present epoch is identified as 
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Figure 6: Constraints on the Hu-Sawicki f{R) gravity model ((SJ from the RSD data in the (n, A) plane for three different scales 
of = 10 h-^ Mpc (left), 30 (middle), 60 Mpc (right). The lower-left regions are excluded by the confidence levels 
indicated along with the contours. 



Covariant Galileon Extended Galileon 




Figure 7: Constraints on the covariant Galileon model given 
by the functions ([T2l)-(IT5l) from the RSD data in the {a, P) 
plane for the scale k^^ — 30 /i^^ Mpc. The dashed and dot- 
ted lines represent the regions constrained by the background 
expansion history according to the Union2 and Constitution 
SN la data, respectively. In the grey region (a — 2/3 > 0.505) 
there is an instability associated with negative values of Gcff 
during the transition from the matter era to the accelerated 
epoch. The parameter space that is theoretically valid and 
consistent with background level observations is excluded by 
the RSD data at more than 11 cr CL. 



^DE ~ 0-72 to match with the background constraints 
derived in Ref . [6ll . 

In Fig. Owe plot the observational bounds in the (a, (3) 
plane for 30 Mpc. For the scales ranging 10 Mpc 



Figure 8: Constraints on the extended Galileon model given 
by the functions (l25l)-(l28l) from the RSD data in the (a, /3) 
plane for the scale k^^ — 30 Mpc. The grey scale region 
is theoretically excluded (see Fig. |4|. The lower- left region 
separated by the 95.4 % contour is excluded at more than 
95.4 % CL, and other contours have similar meanings. 



< < 60 Mpc we confirm that the constraints 

are similar to those shown in Fig. [5] Unlike the covariant 
Galileon there is a parameter space in which the model 
is within the 2a observational contour. 

From Fig. |8] we find that the models with larger val- 
ues of a and /3 tend to be favored due to the suppressed 
growth of matter perturbations in high redshifts. The 
case (iii) shown in Fig. [5] corresponds to such an exam- 
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pie, in which case the model is inside the 90.0 % CL 
observational contour. The case (i) in Fig. [5] is excluded 
at 99.0 % CL because of large values of /mcrg in both 
high and low redshifts. The case (ii) is outside the 95.4 
% CL contour due to the large values of fmas in high 
redshifts in spite of the suppressed evolution of fmas for 
z < 0.5. 

It should be noted again that the observational bounds 
in Fig. [8] correspond to the model parameter s — 0.2. 
For the s values close to the growth rate of matter 
perturbations gets smaller and approaches the ACDM 
rate. Obviously such models should also be compatible 
with the RSD data. Note that the cross-correlation be- 
tween the large-scale structure (LSS) and the integrated- 
Sachs- Wolfe (ISW) effect in CMB can also provide addi- 
tional constraints on the extended Galileon model. The 
joint data analysis of the RSD and the LSS-ISW cross- 
correlation is beyond the scope of our paper. 

V. CONCLUSIONS 

In this paper we have placed observational constraints 
on dark energy models based on f{R) gravity, covari- 
ant Galileon, and extended Galileon from the latest data 
of galaxy rcdshift surveys (WiggleZ, SDSS LRG, BOSS, 
and 6dFGRS). In these models the General Relativistic 
behavior can be recovered at short distances under the 
chameleon mechanism or the Vainshtein mechanism. On 
scales relevant to large-scale structures the modification 
of gravity manifests itself for the growth rate of matter 
density perturbations. This growth rate is related to the 
peculiar velocities of galaxies, which can be constrained 
from the redshift-space distortions. 

As an explicit example we studied the f{R) model (O 
proposed by Hu and Sawicki, which is general enough to 
cover basic properties of f{R) dark energy models. At 
the level of perturbations this model mimics the ACDM 
for the rcdshift z larger than the critical value Zc given 
by Eq. (|50p . but the deviation appears for z < Zc- For 
smaller values of n and A the transition redshift z^ gets 
larger, so that the growth of matter perturbations is more 
significant at late times. 

We placed observational bounds on the Hu-Sawicki 
f{R) model in the (n. A) plane. Since the growth rate in 
this model depends on the scale of density perturbations, 
we derive constraints assuming three different wavenum- 



bers of = 10, 30, and 60 Mpc which are relevant 
to the galaxy redshift surveys considered here. As we see 
in Fig. ini the constraints on model parameters tend to 
be tighter for smaller scales. For = 60 Mpc the 
models with n < 2 and A < 3 are outside the 2a observa- 
tional contour, whereas for k^^ = 10 Mpc the models 
with 71 < 2 and A < 8 are excluded at the 2a CL. These 
results show that the recent RSD data do not favor the 
large deviation of fmas from that in the ACDM, whose 
property can be confirmed in Figs. [1] and [H 

For covariant Galileon there is a tracker solution along 
which the equation of state wbe of dark energy changes 
from —2 (matter era) to —1 (de Sitter era). The likeli- 
hood analysis of SN la, CMB, and BAO shows that only 
the late-time tracking solution is allowed from the data. 
Using the parameter space constrained from the back- 
ground cosmology, we found that the covariant Galileon 
model is excluded at more than lOcr CL. This is asso- 
ciated with the fact that the cosmic growth rate for the 
covariant Galileon is much larger than that in the ACDM, 
as confirmed in Fig.|31 In the range 10 Mpc < k~^ < 
60 h^^ Mpc the evolution of fm^s is practically indepen- 
dent of the scales. 

In the extended Galileon scenario the equation of state 
for the tracker is given by wde = — 1 — s during the 
matter era, where s is a positive constant. For the model 
with s = 0.2 we found that there is a parameter space 
which is inside the 2a observational contour. As we see 
in Fig. [HI the models with larger values of a and /3 tend 
to be favored from the RSD data. 

We have thus shown that the recent observations of 
redshift-space distortions are already quite powerful to 
constrain a number of modified gravity models. In up- 
coming future we should obtain more fas measurements 
that are more precise and/or at different redshifts by 
larger data of ongoing surveys or near-future projects 
such as Subaru/FMOS, Subaru/PFS, HETDEX and so 
on. Wc hope that we will be able to approach the origin 
of dark energy in foreseeable future. 
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